Photon- Induced Magnetization Changes in Single-Molecule Magnets 
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Microwave radiation applied to single-molecule magnets can induce large magnetization changes when the 
radiation is resonant with transitions between spin levels. These changes are interpreted as due to resonant 
heating of the sample by the microwaves. Pulsed-radiation studies show that the magnetization continues to de- 
crease after the radiation has been turned off with a rate that is consistent with the spin's characteristic relaxation 
rate. The measured rate increases with pulse duration and microwave power, indicating that greater absorbed 
radiation energy results in a higher sample temperature. We also performed numerical simulations that quali- 
tatively reproduce many of the experimental results. Our results indicate that experiments aimed at measuring 
the magnetization dynamics between two levels resonant with the radiation must be done much faster than the 
>20- ps time scales probed in these experiments. 

PACS numbers: 75.50.Xx, 61.46.-l-w, 67.57.Lm, 75.45.-(-j 



Single-molecule magnets (SMMs) have been intensively 
studied for more than a decade. There are multiple moti- 
vations for studying these systems. SMMs are some of the 
smallest bistable magnetic units, exhibiting hysteresis below 
a blocking of a few Kelvin. With each molecule in a crystal 
nominally chemically identical, they are better characterizable 
than almost all other systems of nanomagnets. From a prac- 
tical point of view, SMMs represent the ultimate in magnetic 
storage density. At the same time, SMMs exhibit unique quan- 
tum signatures, such as tunneling between spin-orientation 
states (i.e. m levels)'-^-^ and interference between tunneling 
pathsii 

In the last few years, several groups have been studying the 
interaction of SMMs with radiation in an attempt to control 
the magnetization dynamics^'^iiSiSiifliiiii^. Beyond electron- 
spin resonance (ESR) on SMMs,- '- these experiments mea- 
sure changes in the sample magnetization in the presence of 
continuous and pulsed microwave radiation. Ideally, the use 
of radiation to manipulate the dynamics has the advantage that 
population is moved between only two levels - those resonant 
with the applied radiation - while other control parameters, 
such as temperature, affect the population dynamics of multi- 
ple levels. Motivation for these experiments comes from the 
possibility that these systems might be used for quantum com- 
puting, as has been explored theoretically.'^ For that purpose, 
it is essential to measure the lifetime Ti of excited states, as 
well as their decoherence times T2, both of which can in prin- 
ciple be measured with time-domain techniques. Measuring 
these quantities will allow for a deeper understanding of how 
the molecular magnets interact with their environment. For 
example, spin-phonon interactions most likely limit the life- 
time of the excited states. However, these interactions are still 
not fully understood. Calculations of the relaxation rate of 
the molecular magnets using standard spin-phonon interac- 
tions yield results that are orders of magnitude slower than 
measured.'* On the other hand, including interactions that al- 
low a phonon to induce Am = 2 changes in spin give results 
that are in reasonable agreement with the data;'^ '^ however. 



the existence and amplitude of such interactions have been de- 
bated and other mechanisms have been proposed.''-^" Thus, 
a direct measurement of the excited-state lifetimes may help 
shed light on the actual spin-phonon relaxation mechanisms in 
the SMMs. More generally, an understanding of decoherence 
in the SMMs will help address the question of how objects 
lose their quantum properties in the macroscopic limit. 

In this paper, we present the results of studies of the mag- 
netization dynamics of SMMs in the presence of resonant mi- 
crowave radiation. Our results indicate that large magnetiza- 
tion changes can be induced in SMMs by the application of 
resonant microwaves. However, much of these changes are 
due to resonant sample heating by the microwaves on the time 
scale of the spin's characteristic relaxation time. 

The spin dynamics of SMMs can generally be described by 
the effective Hamiltonian 



H = -DSl - BSt - gfiBS.H, + H' 



(1) 



where H' contains terms that do not commute with Sz and 
therefore produce tunneling, and D and B are anisotropy con- 
stants. The first term impels the spin to align with or opposite 
the (easy) z axis, giving rise to a double-well potential, illus- 
trated in Fig. [2 The energy levels are then approximately the 
eigenstates of Sz, with magnetic quantum numbers m, as in- 
dicated in the figure. When a magnetic field Hz is applied 
along the easy axis, it tilts the potential, lowering the energy 
of, say, the right "up" well while raising the energy of the op- 
posite well. Resonant tunneling occurs when the field causes 
levels in opposite wells to align. 

Many of our studies have focused on the so-called Feg 
molecular magnet (Fe802(OH)i2(tacn)6) for which S = 10, 
D = 0.292 K and B is negligible. When continuous wave 
microwave radiation is applied dips occur in the magnitude of 
the equilibrium magnetization whenever the field brings a pair 
of neighboring levels into resonance with the radiation. These 
results represent a variation on traditional ESR spectroscopy. 
In our first experiments, using low-amplitude radiation, we 




FIG. 1: Double- well potential for a typical SMM. The right well 
corresponds to the spin pointing "up" and the left corresponds to it 
pointing "down". A magnetic field tilts the potential, making one 
well higher than the other. The energy levels shown are for the 
Feg SMM but the level distribution is qualitatively similar for other 
SMMs. Mechanisms by which magnetization can be changed are il- 
lustrated by the arrows. Tunneling is generally only significant when 
levels in opposite wells are aligned, as in the figure. 



found very small changes (less than 1%) in the magnetiza- 
tion for the lO-to-9 transitioni^ In the current study, we used 
a high-power backward-wave oscillator and a resonant cavity 
consisting of a two-inch piece of WR-6 waveguide terminated 
on both ends with gold or copper foil. A small coupling hole 
at one end allowed radiation to couple into the cavity from 
gold-plated stainless steel WR-10 waveguide that piped the 
microwaves into the cryostat. Depending on the positioning 
of the cavity caps and the mode excited, the Q could be as 
high 2500. The sample was mounted on a thinned side wall 
of the cavity, opposite a Hall sensor located outside the cavity, 
and oriented with its a axis parallel to the field direction, re- 
sulting the easy axis being tilted by ^^16° relative to the field. 
Samples were on the order of 100 /im in the longest dimen- 
sion, much smaller than the wavelengths used. Samples were 
mounted using silicone-based vacuum grease, which has a low 
thermal conductivity. (The sample was cooled indirectly by 
a few Torr of helium exchange gas filling the sample cham- 
ber.) With this new experimental configuration, we observed 
very large photon-induced changes in magnetization;^ an ex- 
ample is shown in Fig. |2j where several transitions are ob- 
served, as indicated. There are two noteworthy features of 
this data. First, the dips are much larger than one would first 
expect: in the limit in which the lO-to-9 transition is saturated, 
at most 50% of the initial ground-state (m =10) population is 
transfeiTed into the excited (m = 9) state, producing a max- 
imum 5% decrease in the magnetization, while the observed 
dips are much larger Second, we observe dips that correspond 
to transitions between excited states. The sample temperature 
when the microwaves are not on resonance with a transition is 
^4K, smaller than the energies of the excited states. Thus, the 
populations of the excited states should be small and all of the 
transitions between them should have a significantly smaller 
amplitude than the lO-to-9 transition. Instead, all of the dips 
have similar amplitudes and all are larger than would be ex- 
pected even if the transitions were saturated. 

We surmise, therefore, that magnetization changes must in- 
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FIG. 2: Magnetization as a function of field for the Feg SMM both in 
the presence of 121.67-GHz cw microwave radiation and with no ra- 
diation. The radiation heats the sample, especially when it is resonant 
with a transition between spin energy levels, producing the observed 
dips. The cavity Q is ~2000 and the power incident on the cavity 
(after accounting for losses in the waveguide and other components) 
is -^3 mW. 



volve population being moved to levels other than the two 
involved in the radiative transition. In fact, it is likely that 
population is moving to the opposite metastable well. The 
simplest explanation for our observations is that the radiation 
heats the sample significantly when it is on resonance with a 
transition between levels. The dips in magnetization then rep- 
resent the simple fact that the equilibrium magnetization de- 
creases at higher temperatures. We have seen direct evidence 
for such heating using a thermometer to monitor the cavity 
temperature and observed an increase in temperature when- 
ever one of the magnetization dips occurs.'' The fact that the 
magnetization dips occur even for transitions between excited 
states we believe is due to a bootstrapping process in which 
radiative transitions heat the sample, increasing the excited- 
state populations and, in turn, allowing for more photons to 
be absorbed. 

Further evidence that the observed magnetization dips are 
due to radiative heating is provided by a study of another 
SMM, Mn403Cl(02CCH3)3(dbm)3, or simply Mn4 for short. 
This molecule, with a spin of 9/2 and D = 0.76 K, has a 
smaller barrier than Feg. For these samples, the easy axis was 
closely aligned with the magnetic-field direction. Equilibrium 
magnetization curves taken in the presence of microwaves at 
three different frequencies are shown in Fig. Two sets of 
dips occur One dip moves to higher fields with increasing 
frequency, consistent with a transition between levels in the 
right well, where levels get further apart with increasing field. 
We identify this dip with the 7/2-to-5/2 transition. The other 
dip moves to lower fields with increasing frequency, indicat- 
ing that it belongs to a transition in the left well, which we 
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identify as a transition between m = -9/2 and m = -7/2. Since 
this latter transition occurs in the metastable well, if radiation 
were only transferring population between the two states of 
the transition, the equilibrium magnetization should increase 
in the presence of radiation. The fact that we instead observe a 
decrease in magnetization for this transition indicates that the 
primary effect of the radiation is heating of the sample. 
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FIG. 3: Equilibrium magnetization as a function of field when cw 
microwave radiation is applied to the Mn4 SMM. Data for three dif- 
ferent frequencies is shown, as indicated. The nominal system tem- 
perature is ~ 2.4 - 2.5 K. At each frequency, two dips are observed, 
corresponding to the transitions indicated. 

To better understand how this heating takes place, we per- 
formed a study of the magnetization dynamics when pulses 
of radiation were applied to the Feg SMM. Fig. |3 shows the 
change in magnetization (normalized by the sample's satura- 
tion magnetization) at (nominally) 1.8 K and 0.17 T as a func- 
tion of time during and after 119.77 GHz (corresponding to 
the lO-to-9 transition) radiation pulses of various durations, as 
indicated. The magnetization decreases during the radiation 
pulse, following a "backward-S"-shaped curve, seen for the 
longest pulses. (The data exhibit abrupt jumps whenever the 
radiation is turned on or off, which are due to an instrumental 
artifact.) In addition, except for the longest pulses, the mag- 
netization continues to decrease after the radiation has been 
turned off and eventually reaches a minimum. These effects 
are not observed (or at most are very weak) when the radia- 
tion is not resonant with a transition between spin states.^ We 
interpret the results as being due to an indirect heating of the 
spin system by the radiation. As described previously,^ after 
photons are absorbed by the spin transition, the spins relax 
through the emission of phonons. These phonons then rapidly 
thermalize, heating the lattice to a higher temperature. The 
spin system must then come to thermal equilibrium with the 
lattice, a process that takes the spin's characteristic relaxation 
time for thermal activation and/or thermally assisted tunnel- 



ing. At the temperatures of our experiments, the relaxation 
time for SMMs is known to follow an Arrhenius law 



where U is the effective activation barrier and tq ^ 10^^ s 
for many SMMs including Peg. The backward-S shape can 
be qualitatively interpreted as follows. At short times (< 0.1 
ms), the magnetization does not change very much because 

(a) the amount of energy absorbed by the crystal is small and 

(b) the spins have not yet had enough time to come to thermal 
equilibrium with the new lattice temperature. At somewhat 
longer times ('^0.1 - ~1 ms), the magnetization decreases 
rapidly as the spins try to reach thermal equilibrium with the 
higher-temperature lattice. The spins can relax increasingly 
fast at progressively longer times as the temperature of the 
lattice continues to increase during the pulse. At the longest 
times (> 1 ms), the magnetization begins to saturate as the 
entire system (spins and lattice) reach a well-defined steady- 
state temperature and the rate of radiation energy absorption 
is balanced by the rate heat to be carried away by the cryo- 
stat. For the longest pulse durations, the magnetization does 
not decrease after the radiation has been turned off, consistent 
with the interpretation that the spins are close to equilibrium 
with the lattice at these time scales. 
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FIG. 4: Magnetization change (normalized to the saturation magne- 
tization) as a function of time when pulses of 119.77 GHz radiation 
are applied to the Feg SMM. Data for various pulse durations, as in- 
dicated, are shown. The abrupt jumps when the radiation is turned on 
or off are due to an instrumental artifact. The nominal system tem- 
perature is 1.8 K and the applied magnetic field is 0.17 T, at which 
the lO-to-9 transition is resonant with the radiation. The radiation 
power incident on the cavity was 6.2 mW and the cavity Q was sev- 
eral hundred. 

A more quantitative analysis of the data can be carried out 
using the decrease in magnetization that occurs after the radi- 
ation has been turned off. The inset of Fig.|5]shows a typical 
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FIG. 5: Relaxation rate as a function of pulse duration for the Fes 
SMM. The rates are extracted by fitting the post-pulse decays in 
Fig- El to an exponential function, an example of which is shown in 
the inset. The curve in the main figure is the result of numerical 
calculations, as described in the text. 



such decrease, as well as a fit to an exponential decay. The 
fact that such a fit works well is notable, indicating that dur- 
ing this period, the lattice temperature is reasonably constant 
so that the spins have a single definable relaxation rate. This 
is confirmed by the fact that the time for the magnetization to 
return to equilibrium after the pulse is much longer than the 
timescale of the decay. In previous work, we used such fits to 
extract relaxation rates (l/r) from the post-pulse decays and 
found that the rate is maximum at fields for which the absorp- 
tion for the lO-to-9 transition is strongest or for which reso- 
nant tunneling is active, confirming that the rate corresponds 
to the spin's relaxation across the barrier^ Here we apply a 
similar analysis to the data in Fig. 0]to see how the rate de- 
pends on pulse duration; the results are shown in the main 
part of Fig. |5] There is a clear trend to increasing relaxation 
rate for longer pulse lengths. This is consistent with more 
heating with longer pulses. Extrapolating the data in Fig.|5]to 
zero pulse duration yields a relaxation rate of ~1 kHz. The 
value is roughly an order of magnitude faster than other pub- 
lished measurements of the relaxation rate in this system,^'''^^ 
although these do not report values for our precise experimen- 
tal conditions (1.8 K and 0.17 T). 

Figure |6l shows 0.2-ms-pulsed-radiation data for several 
values of radiation power. The behavior is quaUtatively sim- 
ilar to that in Fig.|3] the magnetization decreases during the 
pulse, then continues to decrease after the pulse, reaching a 
minimum, before slowly returning to equilibrium. At higher 
powers the magnetization changes are larger, but the behavior 
does not change qualitatively. We again analyze the decay af- 
ter the pulse by fitting to an exponential. The relaxation rate 
extracted this way is shown in Fig.0 The rate increases with 



power and has a zero-power intercept of just above 1 kHz, 
similar to what is found from the pulse-duration dependence. 
The observed increase in relaxation rate is in qualitative agree- 
ment with our interpretation: the higher the radiation power, 
the higher the lattice temperature and the faster the spins relax. 
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FIG. 6: Magnetization change as a function of time when 0.2-ms 
pulses of I19.77-GHz radiation are applied to the Feg SMM. From 
top to bottom, the curves correspond to power levels of 0.21 mW, 
0.34 mW, 0.59 mW, 0.76 mW, 1.00 mW, 1.24 mW, 1.48 mW, 1.86 
mW, 2.24 mW, 3.00 mW, 3.96 mW, 4.96 mW and 6.20 mW. The 
cryostat temperature is 1.8 K, magnetic field 0.17 T and cavity Q 
several hundred. 

It is interesting to note that Petukhov et al.}^ also working 
with Feg do not see large decreases in the magnetization after 
the radiation is turned off. The reason for the disparity been 
our results and theirs is unclear, but may have something to do 
with thermal anchoring of the sample. 

We simulated the non-equilibrium dynamics of our system 
by tracking the population of the ground state as it is driven 
out of equilibrium by the heating of the lattice and attempts to 
regain equilibrium through thermal relaxation. We first calcu- 
late the photon-induced transition rate for the lO-to-9 transi- 
tion. This can be done using the standard expression^ 

w = ^\t.\'F{oj), (3) 

where = {guB? |(10| (5+ + 5_)/2 |9)|^ is the square 
of the magnetic dipole matrix element. Hi is the amplitude of 
the radiation field and F{lo) is the lineshape function, which 
has a value of 7r/T2 for the peak of a Lorentzian. We estimate 
Hi inside the cavity to be ~0.06 Oe for an incident power of 
6.2 mW and cavity Q~370. T2 is - 1.7 x IQ-^" s.** (While 
this small value likely represents inhomogeneous broadening, 
the nature of the broadening is immaterial to the present cal- 
culation as long as the transition is not driven near saturation.) 
With these numbers, we calculate a transition rate of w^AQO 
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FIG. 7: Relaxation rate as a function of microwave power extracted 
from tlie post-pulse decays in Fig.|6| The solid curve is the result of 
numerical calculations, as described in the text. 



. Since this value is much smaller than any spin-phonon 
transition rate, we can conclude that photon absorption alone 
does not appreciably change the populations of the states. So, 
we treat the levels within the right well as all being in thermal 
equilibrium with the lattice during the radiation pulse. The 
magnetization dynamics are then driven by the radiative heat- 
ing of the lattice. As the system heats, the (normalized) pop- 
ulation of the ground state pio relaxes towards thermal equi- 
librium through thermal activation over the barrier with a time 
constant t given by Eq.|2] 

dpioM= -(pio-Pig(r))/r, (4) 

where p^q(T) = g-Sio/fesT^^ j-j^g equilibrium population 
of the ground state at temperature T and Z = ^-Ei/ksT 

i 

is the partition function for the spin system. The energy 
eigenvalues Ei are found by numerical diagonalization of the 
Hamiltonian. 

The temperature of the system as a function of time can be 
found through energy considerations. The rate of change of 
the sample's energy can be related to the energy of the photons 
absorbed: 



dW/dt = whwN{pio - pg) - K(r - Tenv) 

= wfuvNpioil - e-(^«-^^")/'=«^) - KiT - Tenv), 

(5) 

where N ^ 1.7 x lO^"* is the total number of spins in the sam- 
ple. The last term represents the heat flow out of the sample to 
the cryostat at temperature Tenv (1.8 K for our experimental 



situation); k is a phenomenological parameter characterizing 
the thermal conductance. The energy gain by the sample can 
be related to its temperature change by 

dW = mCdT dW/dt = mCdT/dt, (6) 

where m 0.5 /Ug is the sample mass. The specific 
heat C has two contributions. The lattice specific heat is — 
Ciat = 2.1 X 10"^T^ J/5A'*. The magnetic specific heat (per 
molecule) is calculated from the usual definition: C,„ag — 
kB-^{T^ ^q!^ ) . C,nag is not strictly defined when the spins 
are out of equilibrium. In reality, at low temperatures only 
states in the lower well are in thermal equilibrium with the 
lattice while states in the other well are inaccessible, effec- 
tively lowering the specific heat until the spin system comes 
into equilibrium. An effective time-dependent specific heat 
has been considered theoretically. ^^■^^•^^•^^ We chose not to in- 
clude such a careful accounting of the magnetic specific heat 
because it both significantly increases the complexity of our 
calculations and results only in a factor of ^2 decrease in the 
specific heat during the time when the spins are out of equi- 
librium. 

Equating Eqs. |5]and|6j results in an equation for dT/dt. 
This equation and Eq.l^can be simultaneously solved numer- 
ically to find pio and T as a function of time. From T{t) we 
can calculate the relaxation rate (l/r) of the systems after a 
pulse of length t using Eq. |2] In doing these calculations, we 
used experimentally determined values of the necessary pa- 
rameters (radiation power, cavity Q, sample mass, etc.). We 
set the energy barrier U to be 22.7 K, the calculated height 
of the classical barrier at a field of 0.17 T. In order to obtain 
a relaxation rate in the absence of radiation of ~1 kHz (the 
intercepts in Figs.|5]and0 we set r = 3 x 10~^ s, an order 
of magnitude smaller than the published valuer?^ We chose 
K = 5/iJ/Ks, which results in good agreement between the 
measured dependence of the relaxation rate on pulse duration 
and the calculated values, as shown in Fig.|5] Given the free- 
dom in these parameters, the agreement is reasonably satisfac- 
tory. Without adjusting these parameters, we were also able 
to simulate the relaxation rate dependence on incident radi- 
ation power and achieve a good agreement with the data, as 
shown by the solid curve in Fig.0 Having fixed the param- 
eters, we calculated pio and T as a function of t. The results 
of these calculations are shown in Fig.|S] These results show 
that the (lattice) temperature rises quickly and then saturates at 
^2.5 K. More interestingly, pio shows qualitative agreement 
with the data in Fig.0] exhibiting the characteristic backward- 
S shape and reaching saturation on a time scale of ~1 ms. The 
maximum relative change in pio is 12.5 %, which is roughly 
half the maximum relative change in magnetization in Fig.0] 
as expected since AAI/Msat — ^{pio — P-10) — 2Apio. 
For reference, the figure also shows p^^ as a dashed curve. 
The difference between pio and p^Q indicates that the spin sys- 
tem is out of equilibrium for ~ 1 ms, again in fair agreement 
with with our observations. A more quantitative comparison 
between the data and calculations would require inclusion of 
a proper treatment of the specific heat in the calculations as 
well as a better determination of some relevant experimental 
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parameters, such as the relaxation rate in the absence of radi- 
ation. 
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FIG. 8: Results of numerical simulations for the Feg SMM in the 
presence of radiation resonant with the lO-to-9 transition. The tem- 
perature and ground-state population are shown as a function of time, 
with resonant radiation applied continuously starting at t = 0. The 
dashed curve shows the equilibrium value for the ground state popu- 
lation pIq. 



In conclusion, our results indicate that the magnetization 
dynamics induced by microwaves on the time scales of our 
experiments is dominated by resonant sample heating. The 
timescale for the spin system to respond to such heating coiTe- 
sponds to the spin's characteristic relaxation time. In order to 
observe the pure radiation-induced dynamics between the two 
levels involved in the transition, experiments need to be done 
much faster than this relaxation time. In fact, since phonons 
can only be emitted on the time scale of Ti, all heating ef- 
fects could be avoided by working at time scales shorter than 
this. The Hall sensors used in the work described herein have 
a response time of a few microseconds, which is unlikely to 
be fast enough. An inductive pick-up loop is currently being 
used as a fast detector of magnetization changes. Preliminary 
results of that study are discussed elsewhere in this volume. 
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